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Transformations performing on the dependent and/or the independent variables are an useful 
method used to classify PDE in class of equivalence. In this paper we consider a large class 
of U(l)-invariant nonlinear Schrodinger equations containing complex nonlinearities. The U(l) 
symmetry implies the existence of a continuity equation for the particle density p = |'!^|'^ where the 
current j ^ has, in general, a nonlinear structure. We introduce a nonlinear gauge transformation on 
the dependent variables p and which changes the evolution equation in another one containing 
^ only a real nonlinearity and transforms the particle current j ^ in the standard bilinear form. 

We extend the method to U(l)-invariant coupled nonlinear Schrodinger equations where the 
^ ' most general nonlinearity is taken into account through the sum of an Hermitian matrix and an 

anti-Hermitian matrix. By means of the nonlinear gauge transformation we change the nonlinear 
system in another one containing only a purely Hermitian nonlinearity. Finally, we consider 
nonlinear Schrodinger equations minimally coupled with an Abelian gauge field whose dynamics is 
governed, in the most general fashion, through the Maxwell-Chern-Simons equation. It is shown 
that the nonlinear transformation we are introducing can be applied, in this case, separately to the 
f~i ' gauge field or to the matter field with the same final result. In conclusion, some relevant examples 

CL(| are presented to show the applicability of the method. 

^ • Mathematics Subject Classification (2000): 35Q55, 37K05, 37K35 

^: 

I. INTRODUCTION 

^ ' The Schrodinger equation is one of the most studied topics both from a mathematical and physical point of wiev. 
A particular interest is related to the possible nonlinear extensions of this equation. Just one year after the discovery 
' of the Schrodinger equation, Fermi proposed its first nonlinear generalization pTj . 

, In the following years, many nonlinear extensions of this equation have been proposed in literature in order either 
' to explore fundamental arguments of the quantum mechanics, with the usual linear theory representing only an 
, approximation, or to describe particular phenomenological physical effects. Among the many attempts made to 
generalize in a nonlinear manner the Schrodinger equation we recall the Bialynicki-Birula and Mycielski equation 
[6|, with the nonlinear term —b lii{\tfj\'^) tp; the Guerra and Pusterla model [341] . that with the purpose of preserving 
the superposition pri ncip le of the quantum mechanics introduced the nonlinear term (Aj-;/'!/!?/;!) V'; more recently, the 
Weinberg model |95l . |96| , with the introduction of homogeneous nonlinear terms in order to save partially the same 
fundamental principle. 

^ On a phenomenological basis we recall the well known cubic Schrodinger equation [s^ [ssl TSj , used in the study of 
- - - the dynamical evolution of a Boson gas with a (5-function pair-wise repulsion or attraction 4!| and in the description 
of the Bose-Einstein condensation of alcali atoms like ''Li, ^"^Na and ^^Rb 89, 9l]; the model introduced by Kostin 
57, 58, 86, 87], with the nonlinear term i \n{4>/ip*)ip + i {ln{tp/ip*)) tp used to describe dissipative systems, and others 
Ill^l31iiM,i69i]- 

Many nonlinear Schrodinger equations (NLSEs) contain complex nonlinearities. For instance, a nonlinearity of the type 



ai|i/'| V + 02|i/'| V' + * a3^(MJ 4') + ili + i 0.5) dxl'tpr ^ was introduced to describe a single mode wave propagating 



in a Kerr dielectric guide |23|,|2^; the nonlinearity ai Itp]"^ ip + i a2 tp + i Oxxi^ + icn IV'P V'j proposed in i66|,i67i] to take 
into account of pumping and dumping effects of the nonlinear media, is used to describe dynamical modes in plasma 
physics, hydrodynamics, and also solitons in optical fibers ('65'] and references therein); the nonlinearity ai JV'l^V' + 
i 0'2 dxxx'4' + * CI3 ip)-\-iai ^xJV'Ifj^ introduced to describe the propagation of high power optical pulses in ultra- 

short soliton communication systems pa. |43 - I47l [61I [88| . In [l^llJj a NLSE with the complex nonlinearity k {ip* ^ip — 
i/j 'Vip*) 'V ^p + {k/2) 'V{ip* 'Vip — ip V'0*) ip has been introduced to take into account a generalized Pauli exclusion- 
inclusion principle between the quantum particles constituting the system. In (48j . in the stochastic quantization 
framework, starting from the most general kinetic containing a nonlinear drift term and compatible with a linear 
diffusion term, a class of NLSEs with a complex nonlinearity was derived whilst, recently [83 - [84| . a wide class of 
NLSEs has been obtained starting from the quantization of a classical many body system whose underlying kinetic is 
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described by a nonlinear Fokker-Planck equation associated to a generalized trace-like entropy. Finally, the Doebner- 
Goldin equation [iM US, 113 introduced from topological considerations as the most general class of NLSEs 
compatible with the linear Fokker-Planck equation for the probability density p = IV'P, where the nonlinear term 
was derived from the unitary group representation of the infinite-dimensional diffeomorphism group proposed as a 
universal quantum kinematical group [30|. 

In the recent years, an increasing interest has been also addressed to systems of coupled nonlinear Schrodinger 
equations (CNLSEs) , particularly after the invention of high- intensity lasers which have allowed [71*1 the experimental 
test of the pioneering theoretical works on the optical fibers propagation in long-distance communications [36, 37},. 
In fact, single- mode optical fibers are not really single- mode type since two possible polarizations exist. A rigorous 
study of their propagation requires the use of CNLSEs in order to take into account the evolution of the different 
polarized waves. In 1974 Manakov [68| introduced a CNLSE starting from the cubic NLSE by considering that the 
total field is a superposition of two, left and right polarized fields. When ultrashort pulses are transmitted through 
fibers, CNLSEs with complex high derivative nonlinearities arise [63l. [t^ - It^ . ItqI, [93] . CNLSEs are also employed in 
the study of light propagation through a nonlinear birefringent medium, in systems with nonrelativistic interactions 
among the different kind of particles, in spinor Bo se-E instein condensation, in the description of micro-polar elastic 
solids, among the many [U, i, H, [H, M, IZO, III, M, ■ 

Finally, when coupled with gauge fields, NLSEs are useful in the study of some interesting phenomenologies in 
condensed matter physics. For instance, in the Ginzburg-Landau theory of the superconductivity the cubic NLSE is 
coupled with an Abelian gauge field whose interaction is described by means of the Maxwell equations |55j i)0] . Some 
time, gauge field dynamics can be described by the further Chern-Simons term which confers mass to the field without 
destroy the gauge invariance of the theory. These models have particle-like solutions obeying to a non-conventional 
statistics, named anyons p^. l43l. [osj . which can find an application in the study of high-Tc super-conductors [97, 99[. 
In this work we present, in a unify way, some recent results concerning the gauge equivalence among U(l)-invariant 
nonrelativistic quantum systems described by nonlinear Schrodinger equations containing a complex nonlinearity 
[i^ - is^ , [ssj . We discuss, in particular, some aspects of the gauge transformations of the third kind based on its 
capability to change the NLSE in another one containing only a real nonlinearity (in the case of CNLSE complex 
nonlinearities are replaced by non Hermitian nonlinearities and the gauge transformation changes them in purely 
Hermitian nonlinearities) . It is worthy to observe that nonlinear transformations of the type discussed in the following 
have been introduced previously in literature to relate different families of NLSEs [l^llol, CNLSEs [92, [l^ and, more 
in general, nonlinear PDF [9- 12]. 



In [18| the name of "gauge transformations of third kind" has been coined for the class of the unitary nonlinear 
transformations. They differ from the gauge transformations of first kind, which have constant generators, and those 
of the second kind, which have generators depending on the space coordinate and eventually on time, since, the gauge 
transformations of third kind have generators depending functionally, often in a nonlinear manner, on the fields. 
On the physical ground, they are named gauge transformations because, as stated by Feynmann and Hibbs [22|], in a 
nonrelativistic quantum mechanics, all measurements of observables are always accomplished through a measurement 
of position and time. Thus, quantum theories, for which the corresponding wave- functions give the same probability 
density in space at all time are in principle equivalent [Tsj . In particular, when the wave- functions "0 and cj) are 
related each to the other by a unitary transformation, like in the gauge transformations which we are introducing, the 
densities of probability of position coincide = |0p = p and the fields -0 and (f> describe the same system. In this 
sense, nonlinear gauge transformations permit to classify NLSEs in classes of equivalence. Any member belonging to 
the same class, in spite of its nonlinearity, describe the same physical system. 

After the brief introductory next section about the notations used in this paper we begin, in section 3, by considering 
a wide class of canonical NLSEs with complex nonlinearity in a n -|- 1-dimensional space-time (throughout this paper 
we use unities h = c ^ e ^ 1 and we set m — 1/2) 



describing, in the mean field approximation, the dynamics of a nonrelativistic scalar field ijj conserving the particle 
number N — J pcPx. The real 0] and the imaginary W['0*, ■0] nonlinearities appearing in equation (jl.ip 

are smooth functionals of the fields "0, "0* and of their spatial derivatives of any order. Since the nonlinearity in 
the evolution equation is complex the particle current has, in general, a nonlinear structure which differs from 
the bilinear form of the standard linear quantum mechanics. We introduce the Lagrangian formulation both in 
the wave-function representation and in the hydrodynamic representation and we study the U(l) symmetry, which 
plays a relevant role for the purpose of the present work. Than, we introduce a nonlinear unitary transformation 
■0 — >■ that changes the complex nonlinearity W[ip, + iyV[tp, 0*] in another one iy[0, 0*] which turns out to 
be purely real. As a consequence the new current j ^ assumes the bilinear form of the linear Schrodinger theory. In 
l4ll . [59I l60l . [9^ [93! we can find some examples of nonlinear gauge transformation of the third 





(1.1) 
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kind. All transformations introduced there can be systematically obtained with the method here presented. 
In section 4, we generalize the transformation to CNLSEs by considering the following system of equations 

i ^ + iA* + (^W[i'^, *] +iyV[^'^ * = , (1.2) 

where = {ip-^ , ■ ■ ■ , ^p) is a p-dimensional vector of scalar wave-functions and the nonlinearity A[5'^, ^P] = W^[^t, ^'] + 
«yV[^^, ^] is composed by an Hermitian matrix W — (A + A^)/2 and an anti-Hermitian matrix iW = (A — A^)/2. 
We assume that the system (|1.2p has q < p conserved multiplets which implies the existence of q continuity equations. 
Then, we introduce a nonlinear gauge transformation -> $ which transforms the CNLSE in another one with 
a purely Hermitian nonlinearity A'[^''l', \E'] = Ty[>I'''', '5]. As a consequence, the transformed currents assume the 
standard bilinear form. 

Finally, in section 5, we generalize the method to a class of NLSEs minimally coupled with an Abelian gauge field 
Afj_, where the matter field is described by the following equation 

tDai; + D^tlj+(w[ij*, V-, A]+iyV[i;*, V, A]) V = , (1.3) 

with ip the scalar charged field and D^^ = {Dq, D) denotes the standard covariant derivative. The dynamics of the 
gauge field is provided by the Maxwell-Chern-Simons equation 

^d,F^''+ge-'^'^F.,=j:^ , (1.4) 

where F^j^^ = d^A^~ d^A^^ is the electromagnetic field and j^^ = (p, j'^^) is the covariant current. We stress one more 
that the charged current j has, in general, a nonlinear structure due to the presence, in the evolution equation (jl.Sp . 
of the complex nonlinearity W^[V'*, ii, A\+ i W^['(/'*, i^, A]- The gauge transformation can be applied equivalently to 
the matter field or to the gauge field obtaining the same final result: the nonlinearity in the Schrodinger equation 
(II. SP turns out to be purely real and the expression of the charged current reduces to the standard bilinear form. 
In section 6, we collect some explicit examples to illustrate the applicability of the method whilst, in the conclusive 
section 7, we discuss the possible further development about the nonlinear gauge transformations presented in this 
paper. 



II. PRELIMINARY MATHEMATICAL BACKGROUND 



Let M be a complex n-dimensional smooth manifold labeled by the vector a; = (x^ , . . . , x^). Let F : M — be 
the algebra of the functions on M and F : F ^ IR the algebra of the functionals on F of the type G = J G{x, t) d"x. 
Let ipj{x, t) E F with j = 1, . . . ,p a set of p fields on M, with t a real parameter and we denote by H = , . . . 
a p-dimensional vector on = M x . . . x M. We assume uniform boundary conditions, to guarantee the convergence 
of the integrals, by requiring that all fields and their spatial derivatives vanish quickly on the boundary of A4. 
In the following we consider a nonrelativistic canonical system described by the action A — J C[d] d^x dt, where 
C[fl] € F is the Lagrangian density which depends on the scalar fields € M and their space and time derivatives. 
Hereinafter we use the notation between square brackets G[ip] to indicate the dependence of G on the field V' and its 
spatial derivative of any order. Since the theory is nonrelativistic the Lagrangian contains only time derivatives of 
the first order. 

We introduce the variation of a functional G with respect to fl as 

^G=f#,...,#) , (2.5) 



where the functional derivative can be defined by means of the Euler operator 

5G 



given by 
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with = /{dx^^ . . . dx'-^) and J2[k=o] = Sfc°=o S/^ ■ The sum J2ik ^'^^^ multi-index 7^ = (ii, . . . , i„) with 
< iq < k, J2q iq — k and 1 < <? < n. 

It is easy to show, by using equation (|2.7p . that the Euler operator satisfies the following property 

where B and the components of C belong to F whilst V = (9^ , . . . , 9,J is the n-dimensional gradient operator. 



III. NONLINEAR SCHRODINGER EQUATION FOR SCALAR PARTICLES 

We begin by consider nonrelativistic many body systems of scalar interacting particles described, in the mean field 
approximation, through a very genera family of NLSEs. It is useful to recall the main aspects of the canonical theory 
both in the wave-function formulation and in the hydrodynamic formulation. 

a) Wave-function formulation. 

Let us consider the class of canonical NLSEs described by the Lagrangian density 

jo[r, i'] = U^*^-^^)- 1 wi' - mr, ^] , (3.1) 



2 V 9i 9t , 

which is a functional of the 2-dimensional vector O = (■(/;*, ip). The first two terms are the same encountered in the 
standard linear quantum mechanics whilst the last term is the nonlinear potential describing the interaction among 
the particles of the system. We assume that U[ip*, ip] be a real smooth functional of the fields ip and ip* and their 
spatial derivatives which leave the Lagrangian density (|3.ip invariant under transformations belonging to the U(l) 
group that assures the conservation of the total number of particles. As we will show, this condition imposes a 
constraint on the form of the nonlinear potential U[ip, ip*]. 
By introducing the action of the system 

[ C[^\ ^]d''xdt , (3.2) 
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where the domain of integration is the whole real region TZ — M x ]R, the evolution equation for the vector field ip, 
corresponding to the stationary trajectories of the action, can be obtained from the extremal problem 

M = , (3.3) 

where the variation is performed with respect to the vector ft. According to equation (|2.5p . the Euler-Lagrange 
equations for the fields and ip* are given by 

~ Jr I = (3.4) 

■R 

and its conjugate. We recall that the Lagrangian density is defined modulo a total derivative (null Lagrangian) 
which does not give contribute to the evolution equations because, as stated in equation (|2.8p . the variation of a total 
derivative vanishes. 

After performing the functional derivatives we obtain the following NLSE 

z^ + AV^-f A[^*, V-] =0 , (3.5) 
at 

where A = 9^ -|- . . . -|- 9^ is the n dimensional Laplacian operator and the complex nonlinear term Af?/'* , V'] is given by 

"5i! 



K[iP*, ^p] = -^ / U\ip\ iP]crxdt . (3.6) 
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b) Hydr adynamic formulation. 

In the hydrodynamic formulation we introduce two real fields p and S related to the wave-function through the 
polar decomposition 

■il){x, t) = p^/^(a;, t) exp [i S{x, t)j , (3.7) 

or equivalently 

p{x, t) - \^{x, t)\^ , (3.8) 
^(-.^)-^ogf^) . (3.9) 



2 ° V ^{x, t) 

By defining the action A ^ C[p, S] d"x dt through the Lagrangian density 

C[p, 5] = p - [VSf p - U[p, S] , (3.10) 
at Ap 

where U[p, S] is the nonlinear potential in the hydrodynamic representation, from the variational problem SA — 0, 
where now Jl = (p, S), we obtain two real equations 

— + {VSf + UM~W[p,S]^0, (3.11) 

^ + V-{2pVS) + 2pW[p,S]^Q. (3.12) 
ot 

In equations p. lip and p.l2p Uq[p] = — Ap^/^/p^/-^ denotes the quantum potential 0, Isll whilst the two real func- 
tionals W[p, S] and yV[p, S] are given by 

W[p, S] = --^ / U[p, S]d'^xdt , (3.13) 
op J 
n 

W[p, 5] = 1 C/[p, . (3.14) 



n 



' 4l^^4-.]^ (3.15) 



According to the relation 

5V \5p ' 2p SS 
the quantities VF[p, S*] and yV[p, S"] are related to the nonlinearity A[ip, in 

A[i/., VI = (m/[p, S] + I W[p, 5]) V . (3.16) 

Equation (j3.1ip is an Hamilton- J acohiASks equation for the field S whilst equation (j3.12p describes the evolution 
equation for the field p. The last term in equation p. 121) originates from the nonlinear potential U[p, S\ and is a 
source of particle which, in the general case, destroys the conservation of the number of particles of the system. 
Accounting for equation (|3.16p we can rewrite the NLSE p.Sp in 

i ^ + AV + {W[p, S] + i W[p, 5]) V = , (3.17) 

where the complex nonlinearity A[p, S], expressed in the hydrodynamic fields p and S*, is separated in the real W[p^ S] 
and imaginary yV[p, S] part. 
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A. U(l) Symmetry 



Differently from the linear Schrodinger equation which is U(l)-invariant, the presence of the nonlinearity U[p, S] 
generally breaks this symmetry. In fact, equation (I3.12[) . in general, is not a continuity equation for the field p. In the 
following we study the relevant restrictions to the nonlinear potential U[p, S] so that the Lagrangian p.ip turns out 
to be U(l)-invariant. Such invariance implies, according to the Noether theorem [t^, the conservation of the total 
number of particles by restoring a continuity equation for the field p. 
To begin with, we rewrite equation (|3.12[) in the form 

^ + V-jf^+2pW[p,S]^0, (3.18) 

where 

= (V* VV'- VW*) = 2pV5 , (3.19) 

is the bilinear particles current of the standard quantum mechanics. Making use of equation (I3.14p . equation p.lSp 
become 

§ + V • - ^ 1 U[p, S] d-x dt = 0, (3.20) 
11 



and taking into account the definition of functional derivative ()2.6p we get 



it + 



[fe=i] 

This last equation can be written in the form 
where the nonlinear current j has components 



dp „ . dU 



(j^)^ = 2pd,S+j^ j U[p, S]d^xdt , (3.23) 
because, according to the definition p.7|) . the following relation holds 

= -V- ij^^ j U[p, S]d"xdt] . (3.24) 

Remark that the expression of current p. 231) is always defined modulo the curl of an arbitrary functional G\p, S\ 
which does not give contribute to equation (|3.22p . 

As stated before, the Lagrangian p.ip . for a general nonlinear potential V\p, S\, is not U(l)-invariant. In fact, 
equation p.22p describe a very general kinetics process in which the right hand side plays the role of a source of 
particles. Trivially, if the nonlinear potential U\p^ S\ depends on the field S only through its spatial derivatives the 
right hand side of equation (I3.22p vanishes and it becomes a continuity equation for the field p 



(3.25) 
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The conserved quantity associated to equation (|3.25|) is the total number of particles 

N= [ pd'^x , (3.26) 



M 



where the integral is evaluated on the full real region AI (uniform boundary conditions guarantee the convergence of 
the integral). 

Under the assumption that U[p, S] depends on S only through its derivatives, according to the relation 



S 

SS 



I U[p, S]d-xdt = -V-l I U[p, S] d-x dt j , (3.27) 



we can rewrite the expression of yV[p, 5*] in 



W[p, S] = ^V- J[p, S] , (3.28) 
2p 



{J)^ [P, S\ = — ^ / U\p, S\ d'^xdt , (3.29) 



where ,!T[p, S] is given by 

5id,S) 
n 

and the expression of the current p.23p becomes 

3,^3^ +J[p.S]. (3.30) 

Summing up, the nonlinear potential U[p^ S\ generally breaks the U(l) symmetry of the system p.ip . On the other 
hand, if it depends on the field S only through its spatial derivative, the U(l) symmetry is restored. This can be 
clarified if we taking into account that, under a global U(l) transformation (gauge transformation of first kind) 

^^0 = e"V', (3.31) 

the phase S transforms in 

S^S = S + e, (3.32) 

where e is the constant parameter of the transformation. As a consequence, if S appears in the Lagrangian only 
through its derivatives, the transformation p.3ip does not change the Lagrangian of the system. 



B. Gauge transformation 

We introduce a unitary and nonlinear transformation on the field ij) 

Tp{x, t) -> (/)(a;, t) = U[p, S] il^ix, t) , (3.33) 

whose purpose is to change the NLSE (|3.17|) . which contains a complex nonlinearity W[p, S\ + iyV[p, S"], in another 
one containing only a purely real nonlinearity W[p, S]. As a consequence the current j^, given in equation (I3.30p . is 
transformed in another one j ^ — > j'^"'' having the merely bilinear form of the ordinary quantum mechanics. 
Since the transformation is unitary: lA* — , equation p.33|) does not change the quantity 

p[x, t) = \^{x, i)|2 = \cj){x, i)|2 , (3.34) 

representing the density of probability of position of the system. 
The functional lA[p, S] is defined in 



U[p, S] = exp {i(T[p, S]^ 



(3.35) 
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where the generator cr[p, S] is a real functional which relates the phase S of the field (f> with the phase S of the field ip 

S = S + a[p, S] , (3.36) 

since we define 

0(a;, t) = p^^^ix, t) exp (iS{x, t)^ . (3.37) 

When equation (|3.36p is invertible, we can express the phase 5 as a functional of the fields S and p. 
The expression of the generator cr[p, S] is related to the imaginary part yV[p, S] of the NLSE through 

V<j[p, S] = ^ J[p, S] , (3.38) 
2p 

which defines a[p, S] modulo an arbitrary integration constant. The same equation p.38p imposes a condition on the 
form of the nonlinear potential as it follows from the relation V x V ct = (where V x / means i9. — /. with 
i,j = l,...,n) 

Vxf^MUo. (3.39) 



Equation p. 391) selects the potentials U[p, S] and in this way the nonlinear systems in which we can perform the 
transformation p.33|) . For one-dimensional systems this transformation can be always accomplished. 

By plunging the expression of ip{x, t) = I4^^[p, S] (j){x, t) in the NLSE p.l7p it is easy to verify that it reduces in 
the following evolution equation 

i^+^(t) + W[p, 5] = , (3.40) 
which contains only a real nonlinearity W[p, S] given by 

W[p,S]^W~{Vaf + 2VS-VcT + ^, (3.41) 

where W = W[p, S[p, S]] and a = a[p, S[p, S]]. Because the phase S appears in equation (|3.40|) only trough its 
spatial derivatives, as required by the U(l)-invariance, the arbitrary integration constant arising from the definition 
of (j[p, S\, does not produce any effect and can be posed equal to zero. Finally, the last term in equation (I3.4ip 
can be solved using the Hamilton- Jacobi equation p. Ill) and the continuity equation (|3.12l) in order to reduce the 
nonlinearity W[p, S] in a quantity containing only space derivatives. 

We can easily verify that the continuity equation for the field p, obtained from equation p.40p . is now given by 

l + V-if =0, (3.42) 

where, due to the reality of the nonlinearity yV[p, S], the current 

jf^2pVS, (3.43) 

assumes the standard bilinear form. 

Let us briefly discuss the generalization of this transformation to non canonical systems. Firstly, we observe that for a 
non canonical system the two quantities W[p^ S] and yV[p, S] are not derivable from a potential U[p, S]. In particular, 
the real nonlinearity W[p, S] can assume any arbitrary expression whereas the form of the imaginary nonlinearity 
yV[p, S"], constrained by the continuity equation for the field p, is given by 

mp. S] = ^W- J[p, S] , (3.44) 
2p 

for an arbitrary functional 3^[p, S]. The particle current j ^ is still given through the relation 

j^^2pVS + J[p,S], (3.45) 
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but now the functional J'\p^ S\ is related to yV[p, S\ only through equation (13.441) . Following the same steps described 
for the canonical case, it is easy to verify that the transformation (j3.33p with generator defined in the same way 

Va{p,S\^^J{p,S\, (3.46) 

eliminates the imaginary part of the nonlinearity in the motion equation p. 171) which assumes the expression (I3.40p - 
p.41|) . We observe that, differently from the canonical case, equation (|3.39p now constraints only the form of the 
imaginary part yV[p, S\ whilst the real part 'W\p^ S\ is completely arbitrary. 

Finally, let us remark that when the transformation (j3.33l) is applied to a canonical equation, generally it breaks the 
canonical structure of the theory. Consequently, the new NLSE is no more expressible in the Lagrangian formalism. 
Differently, when the transformation is applied to a non canonical system, the new NLSE can acquire a canonical 
structure. It is not hard to show that this is possible if the transformed nonlinearity is a functional depending 

only on the field p. In fact, the real quantity W is related to a nonlinear potential U through equation ()3.13p and 
because the new nonlinearity is purely real, from equation (j3.14l) it follows that the potential V and consequently W 
cannot depend on the field S. 



IV. COUPLED NONLINEAR SCHRODINGER EQUATIONS 

Physical systems, whose dynamics is described by means of CNLSEs are ubiquitous in nature. They occur, for 
instance, in presence of many interacting particle of different spices or in presence of multi polarized laser beams 
propagating in optical fibers. In this last case, any polarized component of the electric or magnetic field, can be 
considered like a "particle state" since its evolution is describable through a NLS-like equation [36] . 
We observe that for many particle systems of different kind, by denoting with the number of the fcth species, 
many possible combinations of conserved multiplets can be realized. In particular, two relevant physical situations 
are given when: 

a) All the quantities iV^, are separately conserved, which is relevant, for instance, in nonrelativistic systems of multi- 
species where process of transmutation from a species to another one is forbidden. 

b) Only the quantity iVtot = X]fc conserved. Relevant examples are given in the study of light propagation in 
optical fibers where each species describes a polarization mode and only the total intensity of the beam is conserved. 
In the following we introduce a wide class of CNLSEs in the form 

z^+AA* + A[p, S']^' = , (4.1) 

where ^1/ = (?/i^ , . . . , ^/j^), p = (p^ . . . , p^) and S = (S^ , . . . , 5*^) are y dimensional vectors. We denote the operator 

valued matrix by an hat (the lower case letter to[i/| denotes its entries) and use the notation between square 

brackets to indicate the functional dependence on the components of the vector w = (w^ , . . . , ) and on its spatial 
derivatives of any order. Without loss of generality we assume the p x p matrix A in a diagonal form. 
We observe that any system of CNLSEs can be always accommodated in the form given in equation (|4.1|) with a 
diagonal nonlinearity A[p, S]. Such nonlinearity can be separated in an Hermitian matrix W = (A + At)/2 and an 
anti-Hermitian matrix iW — (A — A^)/2. Thus, we can pose A[p, S] — W[p, S] + iyV[p, S], where the diagonal 
matrices W[p, S] and W[yO, S] have purely real entries. Such assumption is only for sake of convenience and does not 
imply any restriction on the form of the nonlinearity. 

We will consider a general situation in which the system (j4.1l) has q conserved multiplets of order , with fc = 1, ...,<? 
and X^fe^'fc — where 1 < 9 < P- In this way, the two particular cases a) and b) quoted previously are recognized for 
q = p and q = 1, respectively. 

Let us organize the fields ipi^ belonging to the vector 'if, in 

*=(V'ii, V'lp^ ; i/-^!, V'2p, ; ) , 

1st multiplct 2nd multiplct gth multiplct (^■'^) 

and, from now on, we relabel the fields ijj. in ^/;^, where, the first index k refers to the fcth multiplct of order p^, , 
whereas the second index I, with I < I < p,., refers to the Ith field inside to the multiplet k. 
Canonical system (j4.ip can be obtain from the Lagrangian density 



(4.3) 
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where the scalar product in the second term is apphed among the gradient operators. The nonhnear potential U[^\ 
is a real functional which depends on the vector fields ^'^ and their spatial derivatives. Accounting for the uniform 
boundary conditions, the potential t/[\l/'l', vanishes together with all its derivatives, at the spatial infinity. 
By introducing the action of the system 



.4 = /£l*.. 



(4.4) 



the evolution equation for the vector field ^I* is given by the stationary trajectories of the action as it follows from the 
variational problem SA = 0, where the variation is performed with respect to the 2p-dimensional vector Jl = (^P^, ^'). 
In this way we obtain the equation 



dt J^-t 



^7[*^ ■f]d"xdt = , 



(4.5) 



n 



and its Hermitian conjugate, which form a system of 2p-nonlinear coupled Schrodinger equations. 
Taking in account of the polar decomposition of the fields ip^^^ in the real fields p^^ and 5^., 



^Jx, t) = pI(^{x, t) exp {iS,,{x, t) 



we can express the variation 5 / Sip*^ as 



In this way, each component of equation (|4.5p can be written in 



U[p, S]d''xdt 



n 



where U[p^ S] is the nonlinear potential in the hydrodynamic representation. Equation 
following matrix form 



z — + A A* 
at 



{W[p, S] + iW[p, S^]) * = , 
where the Hermitian and anti-Hermitian nonlinearities are given by 



(4.6) 



(4.7) 



(4.8) 

can be posed in the 
(4.9) 



W[p, S] = -diag I — / C/[p, S] dxdt 



W[p, S] = -diag 



n 

1 5 



U[p, S] dx dt 



n 



(4.10) 



(4.11) 



Finally, by using the polar decomposition ()4.6p . equation ()4.9p can be separate in a system of 2p nonlinear real coupled 
equations 



^ + ( V5,, r - a,, ^ - w,, [p, S] = 



dPu 
dt 



Pk 

2a,,V-{p,, VSJ + 2p,, wjp, S]={). 



(4.12) 
(4.13) 



The first set of equations (|4.12p is a system of p-coupled Hamilton- Jacobi-like equations for the fields S"^, , whilst the 
second set of equations (|4.13p describes the time evolution of the fields p^., . 
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A. U(l) symmetry 

In the following we consider those systems written in the form (|4.9p which admit a set of q continuity equations 



^ + V-jV,=0, (4.14) 

which assure the conservation of the quantities 

N,^ I dTx . (4.15) 



M 



This imposes some restrictions on the functional dependence of the potential U[p, S] with respect to the fields p and 
S. To obtain such restrictions we recall the following relation 

5 d ^ 5 , . 



so that, by taking in account the expression of the matrix W[p, S], we can rewrite equation (|4.13p in 



dt 

d 

' dS. 



V • [ 2 a,, p,, VS., + J U[p, S] d-x dt 

U[p, S]d''xdt^O . (4.17) 



By summing this equation on the index I, with I < I < p^, we obtain 

dt 



^ + ^- U"l + -^^^ iP^ ^]) + 4 [P, S]^0, (4.18) 



where 



and 



with 



Moreover, we posed 



with 



Pk 
1=1 



il°i=Ee' (4-20) 

1=1 



jI"L =2a,,p,,V^,, . (4.21) 



^.[P, ^l^E-^-t/^^'^] ' (4-22) 

1=1 



iJ.. h [P, S] = j^Q^^ J U[p, S] d-x dt , (4.23) 



whilst 



Pk 



lMS]=Y.IJp,S] , (4.24) 



1=1 
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I mS]d-xdt. (4.25) 



with 

T = - 

By comparing equation (|4.14p with equation (|4.18p we obtain, as a condition, that the functionals [p, S] must be 
expressed as the gradient of a set of functionals [p, S] 

IdpJ] = V-Gdp,S]. (4.26) 

We remark that the expression of the functionals G^. [p, S] is determined univocally from the nonlinear potential 
U[p, S] through equations (|4.24p . ()4.25p and (|4.26p which select, in this way, the class of the Lagrangians (|4.3p of the 
family of CNLSEs compatible with the set of continuity equations (|4.14p . If the conditions (|4.26p are accomplished, 
equations (|4.18p become a system of q continuity equations, where the total currents of the kth multiplet j is given 
by 

j... = + [P' S] + G, [p, S] . (4.27) 

We recall that, as it follows from the Noether theorem, equations (I4.14p are consequence of the invariance of the 
Lagrangian with respect to a global unitary transformation 

<I'^$ = C/^', (4.28) 

where 

f/ = diag(exp(ie)) , (4.29) 

and 




times P2 times times (4.30) 



are the constant parameters of the transformation. 

In fact, the Lagrangian (|4.3p is invariant under the transformation (|4.28p if the nonlinear potential U[p, S] changes 
according to 



6U[pJ]^-j2^.^-GMS], (4.31) 



fc=i 

where G^, [p, S] are arbitrary functionals. We recall that in this way the motion equation (|4.5I) does not change because 
the Lagrangian density is always defined modulo a total derivative of an arbitrary functional. Taking in account for 
the independence of the parameters , from equation (j4.31l) we obtain 

E a|- / ^iP' ^] ^"^ dt^-V-G, [p, S] , (4.32) 

which, according to the definitions (j4.24l) and (|4.25l) coincides with the condition (j4.26p . In addition, because the 
parameters e^. are constants, the potential U[p, S] can depend on the phases S^, also through their spatial derivatives 
of any order. 

B. Gauge transformation 

We are ready to generalize the nonlinear gauge transformation described in the section 3.2 to the family of CNLSEs 
under inspection. Let us introduce the following transformation 

*(a;, t) $(a;, t) = U[p, S] *(a;, t) , (4.33) 
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where U[p, S] is a diagonal and unitary matrix: W =U ^ . This implies 



whilst the phases 5^, are related to fields i^^, through the relation 
since we define 

4>,,{x, t) = p]J^{x, t) cxp (z5(a;, t)) . (4.36) 

Without lost of generality, the matrix lA[p^ S] can be written in 

i2[p, S] = diag exp [i a[p, S]j^ , (4.37) 

where a = (. . . , ct^., , . . .) is a p-dimensional vector with real components. The generators of the transformation 
(7^, [p, S] relate the phase S of the new field $ with the phase S of the old field ^ according to 

S = S + a[p,S]. (4.38) 

We introduce the generators cr^, [p, S] through the relations 

Va,, [p, 5^] = ( J^,, [p, 5] + n,, [p, ^]) , (4.39) 

where T?,^, [p, S*] are arbitrary real functionals related to G^[p^ S], introduced in equation ()4.26p . in 

Y,n,,[p.S\^G,[pJ\. (4.40) 
1=1 



Consistence of equations (|4.39p implies the following constraints 



V X 



— (J^.jp-, 5]+-R,Jp-, 5] 



= . (4.41) 



These equations select the potential J7[p, S] and, through equations (|4.10p - ()4.1ip . the nonlinear system where the 
transformation can be performed. 

We remark that, according to equations (I4.39p . equation (j4.33p defines a wide class of transformations, one for every 
choice of the set of functionals 'K-^^ [p, S]. Each transformation changes the initial system (|4.9p . with the nonlinearity 
W[p, S] +«W[p, 5*], in another one with a purely Hermitian nonlinearity W'[p, S]. 

Preliminarily, we observe that, within the notation (|4.23p and (|4.25p . the matrix yV[p, S\ assumes the expression 

W[p, S] - -diag f [p-, 5] + V • J,, [p, S]) \ . (4.42) 



.2p,, 

In this way, by performing the gauge transformation, equation (14.91) becomes 

I ^ + v4 A$ + {Wt[p, S] + I Wt[p, 5] j $ = , (4.43) 

where 

m [p, S] = diag L,, - a,, ( Va,, )' + 2 a,, V5,, • Va,, + ^) , (4-44) 



and 

.2p, 



>Vt[p,5]=diag(-^J-,Jp',5]) . (4.45) 
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The functionals [p, 5] are given by 

[P, S] = I,, [p, 5] - V • n,, [p, S] , (4.46) 

and fulfill the relations 

g.F,,[p,5]=0 , (4.47) 
as can be verify by using equations (|4.24p . (|4.26p and (|4.40p . 

We remark that, as a consequence of this last relation, equation (|4.43p admits the following set of continuity equations 

^ + V-jtl-0' (4-48) 

where 

Pfc 



•^**.fc ^^•^**,fcz ' (4.49) 



1=1 

and 

jtL =2a„p,, V5,, , (4.50) 

i.e., the nonlinear currents j are transformed in j — !■ t'-'^' which have the standard bilinear form. 
On the other hand, the system (14.431) . with the nonlinearity Wt[p, S] + iWtip, '^l? can be rewritten in 

9$ ^ 

i— + ylA$ + VF'[p, 5]$ , (4.51) 

with a purely Hermitian nonlinearity W' = {W')^ given in the following block-form 

W'[p,S]^diag(wl[p,S]) . (4.52) 
The X matrices W' [p, S] — [p, S\ + C^. [p, S\ have a diagonal part 



Dk[p, S] = diag ( w^, - a^, (Vct.J^ + 2a,, V5,, • Vct, 



9o 



(9i 



(4.53) 

with purely real entries, and an off-diagonal part 

(a) [p,5]=z-^i^i-^e^(^--^-) , (4.54) 

which result to be Hermitian matrices C, — C'^ . 

k k 

We observe that because the Lagrangian (|4.3p is U(l)-invariant, the arbitrary integration constants, deriving from 
the definition (|4.39p . do not produce any effect and can be posed equal to zero. Moreover, the last term in equation 
(|4.53p can be solved using equations (j4.12p - (j4.13p reducing the nonlinearity W'[p, S] in a quantity containing only 
space derivatives. 

The extension of the method to the case of non canonical coupled systems is almost immediate and can be performed 
following the same steps described at the end of section 3.2. 

For a non canonical system the matrix W[p, S] can assume any arbitrary expression whereas the form of the matrix 
W[p, S] is constrained by the existence of the set of the continuity equations for the fields p,,. Without loss of 
generality we can pose 

w., [p, S]^--^ [p, 5] + V • J,, [p, S]) , (4.55) 
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where now the functionals kT^Jp, S] and [p, S] are no more related to the nonUnear potential U[p, S] through 
equations (I4.23P and (|4.25p . The continuity equations (|4.14p require that the functionals I^,[p, S] still fulfill the 
constraints (|4.26p for an arbitrary set of functionals [p, S]. The total currents ^ are given in equation (|4.27p but 

now the functionals [p, S] and JT^, [p, S] are related to the matrix W[p, S] only through equation (|4.55p . Finally, 
we introduce the transformation ()4.33p with generators ()4.39p which eliminates the anti-Hermitian matrix W[p, S] of 
the nonlinearity and transforms the system of CNLSEs in the form given in equation ()4.5ip with only an Hermitian 
matrix W'[p, S] given stiU through equations (|4.52p - (|4.54p . 
Let us now briefly study separately two particular relevant cases: 

a) CNLSEs conserving the number of each species of particles 

We assume p = q with p^^ — I and replace the double index kl ^ k. From equation (j4.5p we obtain the following 
evolution equation for the quantities p^ 



^ + ^-3,.,+h=0, (4.56) 



where the currents j ^ ^ are given by 



j^ ^)^ = 2a, p, d.S, + ^^^^ / U[p, S] d'^xdt , (4.57) 



n 



while the quantities I^. [p, S] assume the expression 



U[p, S]d"'xdt . (4.58) 



Conservation of the single densities p, implies that all the quantities [p. S] must vanish. This requires that the 
potential U[p, S] depends on the phases only through their spatial derivatives as it follows also from arguments 
based on the invariance of Lagrangian under a global unitary transformation. 
Remark that in this case the matrix yV[p, S] assumes the simple expression 

W[pJ]=diag(--^V-JMS]] , (4.59) 



2 P. 

where the functionals •J'^.lp, S] are defined in equations (|4.22l) and (|4.23p . after posing p^ = 1. All the quantities 
[p, S] introduced in equation (|4.26p reduce to constant vectors which, without lost of generality, can be posed 
equal to zero. This implies that all the functionals J-'^lp, S] vanish and after the transformation, the matrix W'[p, S] 
is reduced in a diagonal form given by 



M^'[p,5] =diag 



(4.60) 



which contains now only a purely real nonlinearity since the off-diagonal part C[p, S] vanishes. Remark that, in this 
case, the gauge transformation is univocally defined because the generators are given by VcTj. [p, S] = S ^ [p, S]/2a^, p^. 

b ) CNLSEs conserving the total number of particles 

We pose q— I and replace the double index kl ^ I. From equation (|4.5p we obtain the following evolution equation 
for the density p^^^ 

^ + ^-j + I.JpJ]^0, (4.61) 

where 

p 

1=1 
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is the total density of particles and the current j is given by 



P j S f 

\ n 



U[p, S]d''xdt 



Conservation of p^^^ require that [p, 5] , defined by 



'=1 ' n 



and can be expressed in 

so that equation (j4.6ip becomes a continuity equation 



dt 



V • j = , 

tot ■ 



where the total current j^^^ is given in 

By performing the transformation ()4.33|) with generator ()4.39p . where 

J2^ApJ]^g[pJ] , 

1=1 



(4.63) 



we obtain the new system of CNLSEs ()4.5ip with an Hermitian nonlinearity W'[p, S] 
diagonal part 



D[p, S] = diag 
contains purely real entries whilst the off-diagonal part 

results to be Hermitian. 



0(7 ' 



2 Vft^ 



(4.64) 

(4.65) 
(4.66) 
(4.67) 

(4.68) 

D[p, S] + d[p, S]. The 
(4.69) 

(4.70) 



V. NONLINEAR SCHRODINGER EQUATION COUPLED WITH GAUGE FIELDS 

In this section we generalize the nonlinear transformation to NLSEs coupled with Abelian gauge fields whose 
dynamic is described by means of the standard Maxwell term with the inclusion of the additional Chern-Simons term. 



A. The canonical model 

We consider a class of NLSEs describing, in the mean field approximation, a system of interacting charged particles. 
The model is furnished by the following Lagrangian density 



where the Lagrangian of the matter field Cm is given by 

Lmir, ^, A^] = [r Av- - V' (Av-)*] - mi' - u[p, s, a] , 



(5.1) 



(5.2) 
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with D^, = (9^ + i Afj) the covariant derivative and U[p, S, A] is the nonhnear potential in the hydrodynamic 
representation depending on the abeUan gauge field = {Aq, —A) only through its spatial components. The 
Lagrangian of the gauge field £g assumes the expression 



CM ^~\f,, F^'^ + I s^^'^ A^ F,^ , (5.3) 



where F^^i, = d^A^ — d,yA^ is the electromagnetic tensor, with 9^ = {d/dt, V). 

In the following, Greek indices take the value 0, . . . , n, the Latin indices assume the value 1, . . . , n and denote the 
spatial coordinates. Indices are lowered and uppered depending on the metric tensor r^^^ = diag(l, —1, . . . , —1). The 
Levi-Civita tensor e'^'"', fully antisymmetric, is defined as e°^^ — 1. The parameters 7 and g weight the contribute of 
the Maxwell interaction and the Chern-Simons interaction. We recall that the Chern-Simons term gives contribution 
only when the dynamic of the system is constrained in a manifold with an even number of space dimensions (like in 
the plane) whilst in an odd number of space dimensions it reduces to a total derivative which does not give contribute 
to the motion equation. 
Starting from the action of the system 

J C[i;*, ^, Af,]d''xdt , (5.4) 
n 

the evolution equations for the fields "0, "0* and A^ are obtained by posing SA = where the variation is performed 

with respect to the 3- vector O = (-0, A^). 

The motion equation for the gauge field assumes the expression 

7 9^F^'^+.g£'^^^F.^=j::^ , (5.5) 

where the covariant current = (p, j has spatial components 

(j^^)^ = 2p(a,5 + A,) U[p,S,A]rxdt. (5.6) 

By observing that F'^'^ — — F"^, from equation (|5.5p we immediately obtain the continuity equation for the field p 

f + V-j.,=0, (5.7) 
which assures the conservation of the total charge of the system. 

On the other hand, the evolution equation for the matter field, as it follows from the Lagrangian density (|5.ip . is 
given by 

i Dti^ + r>V + (W[p, S, A]+i W[p, S, A]) V = , (5.8) 
where the real and imaginary parts of the nonlinearity are given, respectively, by 

W[p, S,A]=--^ f U[p, S, A] d'^xdt , (5.9) 
op J 
n 

and 



W[p, 5, A] = - ^ A j u[p, S, A] d"xdt 
n 



(5.10) 



For consistence, equation (|5.8p must admit the same continuity equation (|5.7|) . 

Following standard arguments, by multiply equation (j5.8l) by '0* a-nd taking its imaginary part, we obtain 



dp „ 



2p{yS-A)] -4; / U[p, S, A]d"xdt = , (5.11) 
J dS 



TZ 
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which can be written in 

% + ^-^^.-§sj U[P,S,A]d-xdt, (5.12) 
n 

where the charged current j now becomes 

" •'Ail 

(j^^)^=2p (3.5 + ^0 + ^(1^ Ju[p,S,A]d"xdt . (5.13) 
By comparing this expression with equation (|5.6p it follows that: 

1) the nonlinear potential U[p, S, A] must depend on the field S only through its spatial derivatives so that the right 
hand side of equation (I5.12p vanishes becoming, in this way, a continuity equation for the field p. 

2) the fields V^* and A must be present in the nonlinear potential through the combination V5 — A. In other 
words the Lagrangian of the matter field can be obtained consistently from the Lagrangian of the scalar field p.ip by 
replacing in it the standard derivatives with the covariant ones — ^ = + iA^ (minimal coupling prescription). 
Since, as a required, U[p, S, A] depends only through the quantity VS*— A and its highest spatial derivatives, equation 
(|5.10p can be written in 

W[p, S,A] = ^V- Ja[p, S, A] , (5.14) 
2p 

where the vector ,Ja[p, S, A] is defined in 

{Ja)\p, S, A] = J U[p, 5, A] d'^xdt , (5.15) 

n 

and the charged current (|5.13p assumes the expresion 

j^^=2p{VS-A) + JA[p,S,A]. (5.16) 



B. Gauge transformation 

Firstly, we recall that the system described by the Lagrangian ()5.ip is invariant over a local U(l) transformation 
(gauge transformation of second kind), accomplished both on the fields tp and A^, by means of 

Af,{x, t) Af,{x, t) ~ df,uj{x, t) , 

(5.17) 

'ipix, t) exp {iuj[x, t)^Tp{x, t) , 

where uj{x, t) is a well-behaved function in the sense of e^'^df^di.Lu = 0, with e'"' the anti-symmetric tensor e^" ~ — e''^. 
Remark that, under the transformation (I5.17p . Lagrangian (|5.3p changes according to 

C.^C. + ^e^-'-d.iu^F^r) , (5.18) 

with an extra surface term which does not change the motion of equations for the fields ip and A^. 

Let us now introduce the gauge transformation of third kind as a unitary nonlinear transformation performed only 

on the field tp 

■ip{x, t) <p{x, t) = U[p, S, A] ip^x, t) , (5.19) 

which allows to eliminate the imaginary part yV[p, S, A\ of the nonlinearity in the evolution equation (15.81) and reduces 
the charged current to the standard bilinear form 

3^^ [p, 5, A] ^ [p, S,A]^2p {VS - A) . (5.20) 
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The unitary functional h([p, S, A] is given by 

Kip, S, A] = exp (i a [p, 5, A] ) , (5.21) 
where the real generator of the transformation a [p, S, A] defined according to 



V<j[p,S,A]^^Ja[p,S,A], (5.22) 
z p 



^^, Ja[p, S, -4] , ^j._23) 



are constrained by 



P 

By performing the transformation (I5.19p . from equation (|5.8p we obtain the following NLSE for the charged field 4> 

iDt(j) + D^(t) + W[p,S,A]^ = 0, (5.24) 
where the real nonlinearity W[p, S*, A] assumes the expression 

W[p, S,A]=W- (Va)^ + 2 (VS - A) -Va + ^ , (5.25) 

ot 

with W = W[p, S[p, S, A], A] and cr = a[p, S[p, S, A], A]. The new phase S of the field (p is related to the old 
phase S of the field tp through the relation 

S ^ S + a[p, S, A] , (5.26) 
and because the nonlinearity in equation (j5.24l) is a purely real quantity the continuity equation for the field p becomes 

| + V-i!fi=0, (5.27) 

with the transformed charged current given in equation (|5.20p . 

Since the nonlinear transformation has been accomplished only on the matter field, the evolution equation for the 
gauge field retains formally the same expression given in equation (j5.5p 

jd,F^''+ge-'^^Fr,=j:^, (5.28) 
but with the transformed charged source j'^^ = {p, j^^j)- 

On the other hand, the presence of the gauge field enable us to introduce a transformation on it, leaving the matter 
field unchanged. 

In fact, let us introduce the following transformation 

A{x, t) ^ x{x, t) = A{x, t) - Va[p, S, A] , (5.29) 

where a[p, S, A] is still defined through equation (I5.22p . 
Accounting for Ff^^, — ~F^^, it follows that 

F^,^ EE d^,Ay - d„A^ = - 9^Xm : (5-30) 

whenever a[p, S, A] is a well behaved function fulfilling the relation 

e"" ^^.^,a[p, S,A]=0. (5.31) 

This implies that, for p and v spatial indices, equation (j5.30l) is trivially satisfied as consequence of condition (|5.23l) . 
differently, for p ov v equal to zero, equation (|5.30p implies the following transformation for the component Aq(x^ t) 
of the gauge field 

d 

Aq{x, t) xo{x, t) = Ao{x, t) + ^ cr[p, S, A] . (5.32) 
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By performing the transformation (I5.29[) and (|5.32p in equation (|5.5p we obtain 

^d^F'^^+ge'^^^^Fr^^f^^, (5.33) 

where the new covariant current 7" = (p, \ with 



fulfills the continuity equation 



Differently, from equation (|5.8p it follows 



=2p(V5-x) , (5.34) 



dp ~(o) 



iDtilJ + D'^ilj + W[p, S,x]i^ = ^ , (5.36) 

which has the same form of equation (15.241) but now the covariant derivative is defined in = + ix^, while the 
real nonlinearity becomes 

W[p, S,x]^W^ {Vaf + 2 (V5 - x) • Va + ^ , (5.37) 
WiihW = W[p_,S,A[p,S,x]]s^^'^<J = <y[p.S,A[p,S,x]]- 

In conclusion, it is worthy to observe that if we introduce the nonlinear transformation both on the matter field and 
the gauge filed, by following the prescription given in equation (j5.17p , the evolution equations (j5.5p and (15. 8p are not 
changed in form because the variations due to the matter field are balanced by the variations due to the gauge field. 
Thus, in this case transformation (|5.19p behaviors exactly like a gauge transformation of second kind. 



VI. APPLICATIONS 



To show the applicability of the nonlinear transformation introduced in this paper, we consider some examples for 
the three cases: scalar NLSEs, coupled NLSEs and gauged NLSEs. Some of the examples here discussed are already 
known in literature. We show that the nonlinear transformations introduced by different Authors can be obtained, in 
a unified way, with the method presented in this work. 



A. Scalar NLSEs 



Let us consider, as a first example, the following 1-dimensional NLSE 

dip d^ip 
dt dx^ 



+ zaji^l^^^ + ia^^ij^ ^0 , (6.1) 

where a^, a^, and a_, are real constants. After introducing the hydrodynamic fields p and S we can write the real 
and imaginary part of the nonlinearity in 

W[p,S]^b,p + b,p^ + b,p— , (6.2) 

and 

= K^, (6.3) 
where — a^, = a^, h,^ = ~ and b^ = (a., + 04)72. 

The canonical subclass of equation (|6.ip is given by posing 63 — —2b^ and admits the following potential 
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Equation (j6.ip conserves the density p and the corresponding particles current is given by 

dS 



J^=2p^ + 6.P^ (6.5) 



After performing the transformation (j3.33p with 

X 

^[p] = ^Jpdx', (6.6) 



equation (|6.ip is changed in 



where 6^ = 6^ - bj2 ~ by 4. 

Equation (j6.1l) contains, as particular cases, some know NLSEs. Among them we recall: 

1) The Chen-Lee-Liu equation [l3| (&i = b^ = 0, ~ —2 6^) which is transformed in the NLSE with real nonlinearity 

W[p,S]=b,p' + b,p — , (6.8) 
ox 

where b^ =36^/4. 

2) The Jackiw-Aglietti equation [1, |4l| (6^ ~ 0, b^ = —3 6^/4 and fe., = — 2 6^) which is transformed in the NLSE with 
real nonlinearity 

W[p,S]^b,p—. (6.9) 

3) The Eckaus equation [l^ (6^ = 63 = 0) which is transformed in the NLSE with real nonlinearity 

W[p] = kp^ , (6.10) 

with 63 — 6.-, — by A. Remark that, when 6^ 7^ we obtain, after transformation, the cubic-quintic NLSE with real 
nonlinearity 

w[p]^b,p+kpy (6.11) 

studied in [13]. 

4) The Kaup-Newell equation ^56,] (6^ — b^ = and b^ — ~3b^/2) which is transformed in the NLSE with real 
nonlinearity 

(99 

W[p,S]^by + b,p— , (6.12) 

with 6, =3 6^/16. 

As a second example we consider the canonical NLSE introduced in [H, HJ] 

dt dx"^ \ dx dx J dx 

where k is a real parameter. The real and imaginary nonlinearities in the hydrodynamic representation are given by 

W[p, S] ^-2kp i 

and 



W[p,S]^-2Kpl — ] , (6.14) 



S] = ^^i(p^'^) . (6.15) 



p dx \ dx 
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They are obtained from the potential 



U[p, 5] = K [p 



dS_ 

dx 



whereas the particles current assumes the expression 



, dS 



By performing the transformation p.33|) with generator 



equation (|6.13p changes in 



5^ 

dt dx"^ 



2kp f dS\ K 9^ log p 
l + np \d^) " 2 ^ dx^ 



= . 



(6.16) 



(6.17) 



(6.18) 



(6.19) 



Remark that although equation ()6.13p can be generalized in any spatial dimensions [5J, [sj] condition p.39p is not 
satisfied in general and the transformation p.33p can be applied consistently only in 1-dimensional case. 
Another example is given by the class of the Doebner-Goldin equations [l^ 



dip 



D 



^ + A^+ ^c.i?Jp, S]+^-R,[p] \ V 



(6.20) 



where the nonlinear functionals R. are given by i?^ = V • {pVS)/ p, = ^pI P- R3 = (VS*)^, R^ ~ VS* • ^p/p and 
i?5 = CVp/p)^. The canonical subclass of equation ()6.20p is obtained for = — ^ D, ~ and = —2 and it 
follows from the potential 



U[p, S]^ DVp-VS + c^^ 



{Vpf 



(6.21) 



The particles current is given by 



j=2pVS + DVp 



and the corresponding continuity equation is the well-known Fokker-Planck equation 

where D is the diffusion coefficient. 

By performing the transformation (I3.33P with generator 



equation (j6.20l) transforms in 



(6.22) 



(6.23) 



3.24) 



i ^ + + ^ 5. i?, [p, 5] = 



(6.25) 



with coefficients — ~ D, c^—c^ — -D/2, C3 = C3, = + [c.^ — 1) D and = — D — (Cj — 1) 

It is easy to verify that the generator (|6.24p satisfies condition (13.39^ and the nonlinear transformation can be 

performed in any n > 1 spatial dimensions. 

As a final example we consider the following family of NLSE 



i^ + Aij+ (W{p, S) + I W(p, S)) 4, 







(6.26) 
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with nonlinearities 



D fj^'^\ 
W{p, S)^--f{p)V-i^\+G[p] 



3.27) 



and 



where 



2p 



/(P) = P 



9 In k{p) 
dp 



and G'[p] is an arbitrary functional of p. Equation (|6.26p can be obtained from the potential 



U[p, S] = -D f{p) Vp-VS + J G[p'] dp' , 



(6.28) 



3.29) 



(6.30) 



and was recently derived in the canonical quantization framework from a classical many body systems described by 

generalized entropies [ssj . 

The particles current is given by 



j=2pVS-Df{p) Vp 



3.31) 



which is the sum of a linear drift current j^^^^ — 2 p V5 and a nonlinear diffusion current j^^f^ = f{p) Vp different 
from Pick's current j'pick = —D'Vp which is recovered by posing k{p) = ap, with a a constant. The diffusive term 
is related to the entropy of the classical system through the relation (with the Boltzmann constant fce = 1) 



S{p) 



cC'x / \nK{p')dp' . 



M 



By performing the transformation p.33|) with generator 



<^[P] 



D 



InK(p) , 



3.32) 



(6.33) 



equation (|6.26|) changes in 



with 



hip)Ap + f2ip) [VpY 



cj) + G[p] = 



3.34) 



hip)^ P y-Q^ InK(p) 
1 a hip) 

which contains a purely real nonlinearity depending only on p. 

In particular, starting from the entropy S = — J^^ p log pd"x, with k{p) — e p, equation (j6.26|) becomes 



P 



2 p 



(6.35) 
(6.36) 



3.37) 



which coincides with the canonical sub-family of the Doebner-Goldin equations described in the previous example. 
After transformation it becomes 



Ap 1 /Vp 

~7~2\y 



, 



(6.38) 
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which was studied previously in (3J|. Remarkably, this equation is equivalents to the following linear Schrodinger 
equation 



(6.39) 



with t: = Vl — where the field x is related to hydrodynamic fields p and S through the relation x = \fp exp(i S . 



B. Coupled NLSEs 



Let us now consider the following 1-dimensional system of CNLSEs 



with nonlinearity 



(6.40) 



dx 



3.41) 



where a. - , /3.^ , 7.^ and e.. are real constants and f- (p) are arbitrary real functionals depending only on the vector 
field p. 

In the hydrodynamic representation the nonlinearity (j6.4ip has Hermitian and anti-Hermitian part given, respectively, 

by 



W[p, S] = diag 
yV[p\ = diag 



Z.ft (^'.1^ + ^.. 9^ 



.1=1 
p 



dx 
P. 9p^ 



fAp) 



El d b. lift -|_ e ^ 
\ " p. dx dx 



(6.42) 
(6.43) 



where 6.. = a. . ~ B. . , c. . — . — e. . , d. . = (a.. + B..)/2 and e. . = (7. . + e..)/2. 

Equation (I6.40[) includes some cases already known in literature. For instance: the vector generalization of the Kaup- 
Newell equation [2^ (a. = 1, c., = 0, —b.. =2d.. = e.. = (3 and f (p) = 0); the coupled Chen-Lie-Liu equation (Type 
I) ^ (a^. = 1, c,^. = e'^ = 0, -V^. = 2dJ = f3, ].{p) = 0); the coupled Chen-Lie-Liu equation (Type II) [92] (a^. = 1, 
h.. = d/= 0, cj'= -2e,^ = /3, /^.(p) ='0); the hybrid CNLSE [11 HI (a. = 1, = 0, -fe,^ = 2d^^ = e,^. = P and 
/j (P) = P J2k Pk)' the vectorial Eckhaus equation [i| (a,^ = 0, (p) =_^k \>fc /°i /'J- Moreover, for p = 2, with 
6 . j -I- 2 = and (p) = / Pi + g p^, f^ip) — 9 Pi + f P2J equation (I6.40p has been studied in [81| . 
The canonical sub-family of equation (I6.40p is given by 6^. — c. — —2d.. = ~2e.. and can be obtained through the 
following nonlinear potential 



dS. 



U[pJ]^-J2KP.P.^ + m 



(6.44) 



where the conditions 5 F{p)/5 p. — f.{p) are assumed. 
We observe that: 

a) when d^. — e.., for i j, equation (|6.40p conserves the densities p. and the currents take the form 



dS. , ,2 ^ 

J*., = 2 a, - (rf„ + ) p^. - 2 ^ d^^ p, p^ , 



p 



(6.45) 



with (p) = -{d^^ + e„. ) p2 - 2 E,^^. d.^ p. p^. and L (p) = 0. 

b) when d.j + e.. — d^. + e.. , equation ()6.40p conserves the total density p^^^ ~ p. , and the total current is given 
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by 



P 

Jtot ^ ^ 



55 



(6.46) 



with [p\ = -{d^^ + e^.^. ) p2 and 



3.47) 



If we choose the functionals TZ-\p\ =Q in the case a) and 



(6.48) 



in the case b), where A.^. = d^. + e.^ , the generators (|4.39l) can be written in the unified form 

1 ^ ? 



(6.49) 



By performing the gauge transformation, from equation (|6.40l) we obtain a new system of CNLSEs for the field $ 
with nonhnearity 



where the diagonal matrix S\ has entries 



(6.50) 



D{p, S\ = diag 



dS, dS 

V. 



P 



i,k=l 



(6.51) 



with 



/i . . = fe . . + A . , 
1/ . = c. . A. 



3.52) 



1 



J 

whereas the off-diagonal matrix C has entries 



(a) [p,S]^^ ^f^ e^(^^-^. 



where 



3.53) 



(6.54) 



We observe that the functionals (|6.54p vanish in the case a) and the nonhnearity W[p, S] reduces to a purely real 
quantity. 

Let us now collect some particular cases belonging to equation (|6.40p . 
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1) Bv choosinff b.. = —X and a. c. =2a.A.. we obtain a system of CNLSEs with a purely real nonlinearity which 
depends only on the fields p. 

^c.^.,p.p,+/^(p)| V^^=0. (6.55) 



When [p) — J2ik ^jik Pi Pk ^ith A^..^ — b.. {b^. — 2 5^.)/4a^. , it reduces to a system of decoupled linear Schrodinger 
equations 

i ^ + a . A(/.^ = . (6.56) 
2) By choosing 6,^. = -A,^ for i ^ j, c,^. = 2 a, A,^. and (p) = X^^fe A^^, with 



= Kk {bkk +3A,,/2) /2 a, , 
K.k^K, +A,j2 + A^.J/2a, , 

^A.,. = K, {K - K,/'^) /2 flfc, for fc 7^ j 7^ i and K 7^ j = i 



3.57) 



we obtain the following system of decoupled Jackiw-like NLSEs 



*^ + «, a^+'?,J.<^.=0> (6-58) 



withr/^ ^{b^^+\^)/2a^. 

3) By choosing b.. = — A^^, Aj.^.^ = c^ - \ .j2a- — A,^ X^jAa^, we obtain the CNLSEs 



dt ' dx^ ^ 

being 77^,, = {c^^ — \.^/a-)/2a^. The nonlinear term in equation (|6.59p has been considered in [loj . 

C. Gauged NLSEs 

Let us consider a system of charged particles undergoing to anomalous diffusion and described by the following 
NLSE 

i Dtip + + A[P, 5, A] V' = , (6.60) 

with nonlinearity 



V • ( V5 -A) Ap ( ^P^ ^ 



+7 — — , (6.61) 
z p 

where = q D, —2a and = a(2g — 3) with a, q and I? constant parameters. Equation ()6.60p must be 
considered jointly with 

7 9;.F'^"^+ffe'^^^F.;.=jL , (6.62) 



describing the dynamics of the gauge field. 

The nonlinearity ()6.6ip can be obtained from the potential 



U[p, S, A] =Z3gp«-i Vp- (VS*- A) + ap29-3lZ^ ^ (6.63) 
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and the charged current j is given by 



with j^°l= 2 p{VS~ A). 

As a consequence the system fulfihs the following continuity equation 

which is a nonlinear Fokker-Planck equation for charged particles. 
By performing the transformation (|5.19p with 

2 q — I 



(6.64) 



3.65) 



(6.66) 



where the integration constant has been chosen to avoid the singularity for q — > 1, equations ()6.60p and (|6.62p are 
transformed in 



Ap 
P 



with /3 = 2 a - g2 1)2^2 and 



where j-^ = (p, j(fj) with = 2 p (V5 - A). 

Similar equations can be obtained equivalently by means of the transformation 

X = A-^p^-'Vp, 



Dq 
2 



(6.67) 
(6.68) 

(6.69) 
(6.70) 



It is worthy to observe that equation (|6.60l) . in the q ^ 1 limit, reduces to the gauged canonical subclass of the 
Doebner-Goldin family discussed in section 6.1 



DV-{VS-A) + 2a^-a(^ 
P \ P 

D Ap 
2 p 



which is obtainable from the potential 



U[p, S, A] = DVp-{V S - A) 



iVpf 



(6.71) 



3.72) 



and the continuity equation (j6.65l) reduces to the linear Fokker-Planck equation for charged particles 



dp 
dt 



V-j^°]+DAp = 



In the same limit the gauge transformation has generator 



and reduces equation (|6.7ip to 



iDt(l) + D^(f> + (3 



P ^\pj 



= , 



3.73) 



3.74) 



3.75) 



with 13 = 20-0^/2. 
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VII. CONCLUSIONS AND COMMENTS 



In this paper we have considered a class of canonical NLSEs containing complex nonlinearities and describing U(l)- 
invariant systems. We have introduced a unitary and nonlinear transformation ip ^ (/) which reduces the complex 
nonlinearity in a real one and at the same time transforms the quantum particles current in the standard bilinear 
form. We have extended the method to U(l)-invariant CNLSEs. For these systems we have generalized the gauge 
transformation with the purpose to change the initial nonlinearity in another one purely Hermitian. It has been shown 
that there are many different possibilities to define the generator of the transformation. For any choice we obtain a 
new CNLSE with a different, but Hermitian, nonlinearity. Finally, we have specialized the method for NLSEs mini- 
mally coupled with an Abelian gauged field. We have shown that there are two different ways to reduce the complex 
nonlinearity in a purely real one: or by a nonlinear unitary transformation on the matter field or, alternatively, by a 
nonlinear transformation on the gauge field. 

In the following let us make some considerations about the transformation studied in the present work. 
Firstly, the problem of the integrability of a nonlinear evolution equation is one of the most studied topics in mathe- 
matical physics. Let us consider the most general U(I)-invariant scalar NLSE in the hydrodynamic representation 

^ + V-{2pVS + J)^0 , (7.1) 
— + {VSy + U^-W = 0. (7.2) 



In the Calogero picture |8l. llll - ll3l |. the system of equations (|7.ip and (|7.2|) is C-integrable if there exists a transformation 
of the dependent and/or independent variables: t T, x X, p R, S S which changes equations (17. ip . (|7.2p 
in 

— + V.(2i?V5)=0, (7.3) 

f)C 

— + {VSr + U,^0, (7.4) 

where V and Uq are the gradient and the quantum potential in the new variables. Equations (|7.3p and (|7.4p constitute 
the well known hydrodynamic representation of the standard linear Schrodinger equation. 

On the other hand, the transformation on the field — > 5 introduced in this paper, reduces the continuity equation 
(|7.ip in the standard form given by equation (|7.3p and can be seen as a first step in the Calogero program. 
Secondly, the most general gauge transformation of the kind discussed in the present work can be stated as 

Tp{t, x) <p{t, x) ^ U[p, S] 'ipit, x) , (7.5) 

which is an infinite dimensional unitary representation of the diffeomorphism group with 

U[p, S]^exp(iuj[p, S]). (7.6) 

As matter of fact the real generator uj[p, S] could be any arbitrary functional depending on the fields p and S. 
For instance, in (isj the generator of the transformation has been assumed in 

S) = ^\ogp+ (A(t) -1})S + e(t, x) , (7.7) 

which produces a group of transformations mapping the Doebner-Goldin equation in itself. We observe that the one 
parameter subclass of this transformation with 9{t, x) = 0, X{t) = 1 and j{t) = constant^ coincides with those studied 
in this work. 

Throughout this paper, the generator of the gauge transformation has been chose with the purpose to make real the 
complex nonlinearity of the NLSE under inspection. Alternatively, nonlinear gauge transformations can be useful 
generalized with the purpose to classify NLSEs in equivalence classes. Any equation belonging to the same class, in 
spite of its nonlinearity, is gauge equivalent, by means of equation (17. 6p . to the others equations belonging to the same 
class. 

For instance, let us consider the following family of NLSEs 

i^ + A^ + A[p, 5] V = , (7.8) 
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with complex nonlinearity 



A[p,S] = f,ip)AS + fJp)Vp-VS + f,{p){Vp)^ 

+ /,(p)Ap+^v(/,(p)Vp) , (7.9) 

where the expression of the imaginary part guarantees the existence of a continuity equation for p. 
The quantities /. (p) are functional parameters fixing the NLSE. Any NLSE belonging to the family of equations (|7.8p 
can be determined univocally through the vector / = {/^(p), . . . , /^(p)}. By performing a gauge transformation with 
generator w(p) depending only on the field p, equation (|7.8p changes in 



with 



z^ + A0 + A[p, 5]0 = O, 



A[p, S] = A (p) AS + h (p) Vp • V5 + A (p) ( Vp)^ 



/,(p)Ap + -V(/,(p)Vp 



(7.10) 

(7.11) 
(7.12) 



What is important is to note that the transformation maintains the same structure in the nonlinearity through the 
presence of the functional groups AS", VS* • Vp, (Vp)^ and Ap, whilst the expressions of the new parameters /. (p) 
are given by 

A=A-2p^, 



dp 



dp 



+ A-2p 



duj\ d 



dp J dp 



doj 
'dp 



duj 
dp 



By eliminating w(p) among these equations we obtain a set of gauge invariants relations 



(7.13) 



fs - fs 



+ f — 
dp ^ dp 



/1 + 2/J /, 



(7.14) 



Given two NLSEs belonging to the family (|7.8p . labeled by the respective vectors / = {/^(p), ■ . ■ , f^ip)} and /' = 
{flip), ■ • . , f^ip)}, if the functionals /^(p) and /^(p) fulfil the relations (I7.14p . the two NLSEs are gauge equivalents 
since there exist a generator uj{p) such that, by means of equation (|7.5I) . transforms the first NLSE, labeled by the 
vector /, in the second NLSE labeled by the vector /'. 

In particular, observing that the linear Schrodinger equation is represented by the vector / = {0, 0, 0, 0, 0}, it follows 
that any NLSE fulfilling the relations 

/. = , 



P 



dp 



(7.15) 
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is gauge equivalent to the linear Schrodinger equation. This sub-family can be linearizable by means of the transfor- 
mation (|7.5p with generator uj{p) — J^{f^{p')/p') dp' . In this sense, equations (|7.15p define the subclass of the family 
of equations (|7.8|) which are C-integrable. 

In conclusion, we have shown that the transformation introduced in the present work allows us to deal, in a unifying 
scheme, different NLSEs already known in literature, obtaining in a systematic way the transformations introduced 
by various Authors. 

A natural continuation of this work could be performed in several ways: 

1) Extending the method to the case of NLSEs coupled with non Abelian gauge fields, which are relevant, for instance, 
in the study of heavy-quark particle systems. 

2) Extending the method to relativistic nonlinear equations. In this context, in [l6| . a relativistic generalization of 
the transformation introduced in [Tsl has been proposed to generate nonlinear extensions of the Dirac equation. 

3) Extending the method to discrete NLSEs, which are particularly relevant in the study of lattice models in condensed 
matter. 



[1] Ablowitz M.J., Benney D.J.: Evolution of multi-phase modes for nonlinear dispersive waves. Stud. Appl. Math. 49, 
225-238 (1979). 

[2] Aglietti U., Griguolo L., Jackiw R., Pi S.-Y., Seminara D. Anyons and chiral solitons on a line:, Phys. Rev. Lett. 77, 
4406-4409 (1996). 

[3] Agrawal G.P.: Modulation instability induced by cross-phase modulation, Phys. Rev. Lett. 59, 880-883 (1987). 
[4] Barashenkov I., Harin A.: Nonrelativistic Chern-Simons theory for the repulsive Bose gas, Phys. Rev. Lett. 72, 1575-1579 
(1994). 

[5] Berkhoer A.L., Zakharov V.E.: Self excitation of waves with different polarizations in nonlinear media, Zh. Eksp. Teor. 

Fiz. 58, 903-911 (1970); [Sov. Phys. JETP 31, 486-490 (1970)]. 
[6] Bialynicki-Birula I., Mycielski J.: Nonlinear wave mechanics, Ann. Phys. (N.Y.) 100, 62-93 (1976). 

[7] Bohm D.: A suggested interpretation of the quantum theory in terms of "hidden" variables, Phys. Rev. 85, 166-193 
(1951). 

[8] Calogero F., Degasperis A., De Lillo S.: The multicomponent Eckhaus equation, J. Phys. A: Math. Gen. 30, 5805-5814 
(1997). 

[9] Calogero F.; Universal C-integrable nonlinear partial-differential equation in 7i-|- 1 dimensions, J. Math. Phys. 34, 3197- 
3209 (1993). 

[10] Calogero F.: C-integrable nonlinear partial-differential equations in n -I- 1 dimensions, J. Math. Phys. 33, 1257-1271 
(1992). 

[11] Calogero P., Xiaoda J.: C-integrable nonlinear PDES .2, J. Math. Phys. 32, 875-887 (1991). 
[12] Calogero P., Xiaoda J.: C-integrable nonlinear PDES .2, J. Math. Phys. 32, 2703-2717 (1991). 

[13] Calogero P., De Lillo S.: The Eckhaus Pde itpt + ip^:c + 2 V + j-^l* = 0, Inv. Problems 3, 633-681 (1987) (Corri- 

gendum), Inv. Problems 4, 571 (1988). 

[14] Chen H.H., Lee Y.C., Liu C.S: Integrability of non-linear Hamiltonian-systems by inverse scattering method, Phys. Scr. 
20, 490-492 (1979). 

[15] Dodonov V.V., Mizrahi S.S.: Generalized nonlinear Doebner-Goldin Schrodinger equation and the relaxation of quantum- 
systems, Physica A 214, 619-628 (1995). 

[16] Doebner H.-D., Zhdanov R.: Nonlinear Dirac equations and nonlinear gauge transformations, (2003); 
arXiv:quant-ph/0304167 

[17] Doebner H.-D., Goldin G.A., Nettermann P.: Properties of nonlinear Schrodinger equations associated with diffeomor- 
phism group-representations, J. Math. Phys. 40, 49 (1999). 

[18] Doebner H.-D., Goldin G.A.: Introducing nonlinear gauge transformations in a family of nonlinear Schrodinger equations, 
Phys. Rev. A 54, 3764-3771 (1996). 

[19] Doebner H.-D., Goldin G.A.: Properties of nonlinear Schrodinger-equations associated with diffeomorphism group- 
representations, J. Phys. A: Math. Gen. 27, 1771-1780 (1994). 

[20] Doebner H.-D., Goldin G.A.: On a general nonlinear Schrodinger equation admitting diffusion currents, Phys. Lett. A 
162, 397-401 (1992). 

[21] Fermi E.: Rend. Accad. Naz. Lincei 5, 795 (1955). 

[22] Feynmann R.P., Hibbs A.R.: Quantum Mechanics and Path Integrals, McGraw-Hill, New- York, (1965). 
[23] Florjahczyk M., Gagnon L.: Dispersive-type solutions for the Eckhaus equation, Phys. Rev. A 45, 6881-6883 (1992). 
[24] Florjahczyk M., Gagnon L.: Exact-solutions for a higher-order nonlinear Schrodinger equation, Phys. Rev. A 41, 4478- 
4485 (1990). 

[25] Fordy A. P.: Derivative nonlinear Schrodinger equations and hermitian symmetric-spaces, J. Phys. A, Math. Gen. 17, 
1235-1245 (1984). 

[26] Gedalin M., Scott T.C.: Optical solitary waves in the higher order nonlinear Schrodinger equation , Band Y.B., Phys. 
Rev. Lett. 78, 448-451 (1997). 



31 



[27] Ginzburg V., Pitaevskii L.: On the theory of superfluidity, Zh. Eksp. Theor. Fiz. 34, 1240-1245 (1958); [Sov. Phys. JETP 
7, 858-861 (1958)]. 

[28] Gisin L.: Microscopic derivation of a class of non-linear dissipative Schrodinger-like equations, Physica A 111, 364-370 
(1961). 

[29] Goldin G.A.: The diffeomorphism group-approach to nonlinear quantum-systems. Int. J. Mod. Phys. B 6, 1905-1916 
(1992). 

[30] Goldin G.A., Menikofi^ R., Sharp D.H.: Difi^eomorpliism-groups, gauge groups, and quantum-theory, Phys. Rev. Lett. 51, 
2246 2249 (1983). 

[31] Grigorcnko A.N.: Measurement description by means of a nonlinear Schrodinger equation, J. Phys. A: Math. Gen. 28, 
1459-1466 (1995). 

[32] Gross E.P.: Hydrodynamics of a superfluid condensate, J. Math. Phys. 4, 195-207 (1963). 

[33] Gross E.P.: Structure of a quantized vortex in boson systems, Nuovo Cimento 20, 454-477 (1961). 

[34] Guerra F., Pusterla M.: A nonlinear Schrodinger equation and its relativistic generalization from, basic principles, Lett. 

Nuovo Cimento 34, 351-356 (1982). 
[35] Hacinliyan I., Erbay S.: Coupled quintic nonlinear Schrodinger equations in a generalized elastic solid, J. Phys. A: Math. 

Gen. 37, 9387-9401 (2004). 

[36] Hasegawa A., Kodama Y.: Solitons in optical communication, Oxford University Press, New York, (1995). 

[37] Hasegawa A., Tappert F.D.: Transmission of stationary nonlinear optical pulses in dispersive dielectric fibers. I. Anomalous 

dispersion , Appl. Phys. Lett. 23, 142-144 (1973). 
[38] Hisakado M., Wadati M.: Integrable multicomponent hybrid nonlinear Schrodinger-equations, J. Phys. Soc. Jpn. 64, 

408 413 (1995). 

[39] Hisakado M., lizuka T., Wadati M.: Coupled hybrid nonlinear Schrodinger equation and optical solitons, J. Phys. Soc. 

Jpn. 63, 2887-2894 (1994). 
[40] Ho T.-L.: Spinor Bose condensates in optical traps, Phys. Rev. Lett. 81, 742-745 (1998). 
[41] Jackiw R.: A nonrelativistic chiral soliton in one dimension, J. Nonlinear Math. Phys. 4, 261-270 (1997). 
[42] Jackiw R., Pi S.-Y.: Self-dual Chern-Simons solitons. Prog. Theor. Phys. Suppl. 107, 1-40 (1992). 

[43] Jackiw R., Pi S.-Y.: Classical and quantal nonrelativistic Chern-Simons theory, Phys. Rev. D 42, 3500-3513 (1990) 

(Corrigendum), Phys. Rev. D 42, 3929-3929 (1993). 
[44] Karpman V.I., Rasmussen J.J., Shagolov A.G.: Dynamics of solitons and quasisolitons of the cubic third-order nonlinear 

Schrodinger equation, Phys. Rev. E 64, 026614-13 (2001). 
[45] Karpman V.I.:, Shagalov A.G.: Evolution of solitons described by the higher-order nonlinear Schrodinger equation. II. 

Numerical investigation, Phys. Lett. A 254, 319-324 (1999). 
[46] Karpman V.I.: Evolution of solitons described by higher-order nonlinear Schrodinger equations, Phys. Lett. A 244, 

397-400 (1998). 

[47] Karpman V.I.: Radiation by solitons due to higher-order dispersion, Phys. Rev. E 47, 2073-2082 (1993). 
[48] Kaniadakis G., Scarfone A.M.: Nonlinear Schrodinger equations within the Nelson quantization picture, Rep. Math. Phys. 
51, 225-231 (2003). 

[49] Kaniadakis G., Miraldi E., Scarfone A.M.: Cole-Hopf like transformation for a class of coupled nonlinear Schrodinger 

equations, Rep. Math. Phys. 49, 203 209 (2002). 
[50] Kaniadakis G., Scarfone A.M.: Cole-Hopf-like transformation for Schrodinger equations containing complex nonlinearities, 

J. Phys. A: Math. Gen. 35, 1943-1959 (2002). 
[51] KaniadaJiis G., Scaxfone A.M.: Nonlinear transformation for a class of gauged Schrodinger equations with complex 

nonlinearities, Rep. Math. Phys. 48, 115-121 (2001). 
[52] Kaniadakis G., Scarfone A.M.: Nonlinear gauge transformation for a class of Schrodinger equations containing complex 

nonlinearities, Rep. Math. Phys. 46, 113 118 (2000). 
[53] Kaniadakis G., Quarati P., Scarfone A.M.: Soliton-like behavior of a canonical quantum system obeying an exclusion- 
inclusion principle, Physica A 255, 474-482 (1998). 
[54] Kaniadakis G., Quarati P., Scarfone A.M.: Nonlinear canonical quantum system of collectively interacting particles via 

an exclusion-inclusion principle, Phys. Rev. E 58, 5574-5585 (1998). 
[55] Kaper H.G., Takac P.: Ginzburg-Landau dynamics with a time-dependent magnetic field, Nonlinearity 11, 291-305 

(1998). 

[56] Kaup D.J., Newell A.C.: Exact solution for a derivative non-linear Schrodinger equation, J. Math. Phys. 19, 798-801 
(1978). 

[57] Kostin M.D.: Friction and dissipative phenomena in quantum mechanics, J. Stat. Phys. 12, 145-151 (1975). 
[58] Kostin M.D.: On the Schrodinger-Langevin equation, J. Chem. Phys. 57, 3589-3591 (1973). 

[59] Kundu A.: Comments on the Eckhaus PDF i Vt + V'xo: + 2 (|VI^)a: + \i>f = 0, Inv. Problems, 4, 1143-1144 (1988). 

[60] Kundu A.: Landau-Lifshitz and higher-order nonlinear-systems gauge generated from nonlinear Schrodinger type equa- 
tions, J. Math. Phys. 25, 3433-3438 (1984). 

[61] Li Z., Li L., Tian H., Zhou G.: New types of solitary wave solutions for the higher order nonlinear Schrodinger equation, 
Phys. Rev. Lett. 84, 4096-4099 (2000). 

[62] Madelung E.: Quantum theory in hydrodynamical form, Z. Phys. 40, 332-336 (1926). 

[63] Mahalingam A., Porsezian K.: Propagation of dark solitons in a system of coupled higher-order nonlinear Schrodinger 

equations, J. Phys. A: Math. Gen. 35, 3099 3109 (2002). 
[64] Malomed B.A., Stenflo L.: Modulational instabilities and soliton-solutions of a generalized nonlinear Schrodinger equation. 



32 



J. Phys. A: Math. Gen. 24, L1149-1153 (1991). 
[65] Malomcd B.A.: Bound solitons in the nonlinear Schrodinger-Ginzburg-Landau equation, Phys. Rev. A 44, 6954-6957 
(1991). 

[66] Malomed B.A., Nepomnyashchy A.A.: Kinks and solitons in the generalized Ginzburg-Landau equation, Phys. Rev. A 
42, 6009-6014 (1990). 

[67] Malomed B.A.: Evolution of nonsoliton and quasi- classical wavetrains in nonlinear Schrodinger and Korteweg-Devries 

equations with dissipative perturbations, Physica D 29, 155-172 (1987). 
[68] Manakov S.V.: On the theory of two-dimensional stationary self-focusing of electromagnetic waves, Zh. Eksp. Teor. Fiz. 

65, 505-516 (1973); [Sov. Phys. JETP 38, 248-253 (1974)]. 
[69] Martina L., Soliani G., Winternitz P.: Partially invariant solutions of a class of nonlinear Schrodinger equations, J. Phys. 

A: Math. Gen. 25, 4425-4435 (1992). 
[70] Matthews M.R., Anderson B.P., Haljan P.O., Hall D.S., Holland M.J., Williams J.E., Wieman C.E., Cornell E.A.: 

Watching a superfluid untwist itself: Recurrence of rabi oscillations in a Bose-Einstein condensate, Phys. Rev. Lett. 83, 

3358-3361 (1999). 

[71] Mollenauer L.F., Stolen R.H., Gordon J. P.: Experimental-observation of picosecond pulse narrowing and solitons in 

optical fibers, Phys. Rev. Lett. 45, 1095-1098 (1980). 
[72] Nalckeeran K.: Exact dark soliton solutions for a family of N coupled nonlinear SchrOdinger equations in optical fiber 

media, Phys. Rev. E 64, 046611-7 (2001). 

[73] Nakkeeran K.: On the integrability of the extended nonlinear Schrodinger equation and the coupled extended nonlinear 

Schrodinger equations, ,J. Phys. A: Math. Gen. 33, 3947-3949 (2000). 
[74] Nakkeeran K.: Exact soliton solutions for a family of coupled nonlinear Schrodinger equations in optical fiber media, 

Phys. Rev. E 62, 1313-1321 (2000). 
[75] Newboult G.K., Parker D.F., Faulkner T.R.: Coupled nonlinear Schrodinger equations arising in the study of monomode 

step-index optical fibers, J. Math Phys. 30, 930-936 (1989). 
[76] Noether E.: Invariante Variationsprobleme, Nachr. d. Konig. Gesellsch. d. Wiss. zu Gottingen, Math-phys. Klasse, 235 

(1918) (Enghsh translation from Travel M.A.: Transport Theory and Statistical Physics 1(3), 183 (1971). 
[77] Olver P.J.: Applications of Lie Groups to Diff'erential Equations, Springer, New York, (1986). 

[78] Pitaevsku L.P.: Vortex lines in an imperfect Bose gas, Zh. Eksp. Teor. Fiz. 40, 646-651 (1961); [Sov. Phys. JETP 13, 
451-454 (1961)]. 

[79] Radhakrishnan R., Kundu A., Lakshmanan M.: Coupled nonlinear Schrodinger equations with cubic-quintic nonlinearity: 

Integrability and soliton interaction in non-Kerr media, Phys. Rev. E 60, 3314-3323 (1999). 
[80] Ryskin N.M.: Schrodinger bound nonlinear equations for the description of multifrequency wave packages distribution in 

nonlinear medium with dispersion, Zh. Eksp. Teor. Fiz. 106, 1542-1546 (1994); [Sov. Phys. JETP 79, 833-834 (1994)]. 
[81] Sakovich S.Y., Tsuchida T.: Symmetrically coupled higher-order nonlinear Schrodinger equations: singularity analysis 

and integrability, J. Phys. A: Math. Gen. 33, 7217-7226 (2000). 
[82] Scarfone A.M.: Stochastic quantization of an interacting classical particle system, J. Stat. Mech.: Theo. Exp. P03012-I-16 

(2007). 

[83] Scarfone A.M.: Canonical quantization of classical systems with generalized entropies. Rep. Math. Phys. 55, 169-177 
(2005). 

[84] Scarfone A.M.: Canonical quantization of nonlinear many-body systems, Phys. Rev. E 71, 051103-15 (2005). 

[85] Scarfone A.M.:Gauge transformation of the third kind for U(l)-invariant coupled Schrodinger equations, J. Phys. A: 

Math. Gen. 38, 7037-7050 (2005). 
[86] Schuch D.: Nonunitary connection between explicitly time-dependent and nonlinear approaches for the description of 

dissipative quantum systems, Phys. Rev. A 53, 945-940 (1997). 
[87] Schuch D., Chung K.-M., Hartmann H.: Nonlinear Schrodinger-type field equation for the description of dissipative 

systems 3. Frictionally damped free motion as an example for an aperiodic motion, J. Math. Phys. 25, 3086-3092 (1984). 
[88] Shchesnovich V.S., Doktorov E.V.: Perturbation theory for the modified nonlinear Schrodinger solitons, Physica D 129, 

115-129 (1999). 

[89] Shi H., Zheng W.-M.: Bose-Einstein condensation in an atomic gas with attractive interactions , Phys. Rev. A 55, 
2930-2934 (1997). 

[90] Stratopoulos G.N., Tomaras T.N.: Vortex pairs in charged fiuids, Phys. Rev. B 54, 12493-12504 (1996). 
[91] Stringari S.: Collective excitations of a trapped Bose condensed gas, Phys. Rev. Lett. 77, 2360-2363 (1996). 
[92] Tsuchida T., Wadati M.: Complete integrability of derivative nonlinear Schrodinger-type equations, Inv. Problems 15, 
1363-1373 (1999). 

[93] Tsuchida T., Wadati M.: New integrable systems of derivative nonlinear Schrodinger equations with multiple components, 

Phys. Lett. A 257, 53-64 (1999). 
[94] Vinoj M.N., Kuriakose V.C.: Multisoliton solutions and integrability aspects of coupled higher-order nonlinear Schrodinger 

equations, Phys. Rev. E 62, 8719-8725 (2000). 
[95] Weinberg S.: Precision tests of quantum mechanics, Phys. Rev. Lett. 62, 485-488 (1989). 
[96] Weinberg S.: Testing quantum mechanics, Ann. Phys. 194, 336-386 (1989). 

[97] Weinberg S.: Understanding the Fundamental Constitutents of Matter, A. Zichichi (ed.). Plenum, New York, (1978). 
[98] Wilczek F.: Fractional Statistics and Anyon Superconductivity, World Scientific, Singapore, (1990). 
[99] Wilczek F.: Magnetic flux, angular momentum, and statistics, Phys. Rev. Lett. 48, 1144-1146 (1982). 



33 

[100] Yip S.-K.: Internal vortex structure of a trapped spinor Bose-Einstein condensate, Phys. Rev. Lett. 83, 4677-4681 (1999). 



